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Abstract. How can a system in a macroscopically stable state explore energetically more 
favorable states, which are far away from the current equilibrium state? Based on continuum 
mechanical considerations we derive a BoussiNESQ-type equation 

pu = d^a{d^u) ^ Pdlu, a; G (0, 1), /3 > 0, 

which models the dynamics of martensitic phase transformations. Here p > is the mass 
density, pd^il is a regularization term which models the inertial forces of oscillations within 
a representative volume of length ^/]3 and a is a nonmonotone stress-strain relation. The 
solutions of the system, which we refer to as microkinetically regularized wave equation, exhibit 
strong oscillations after times of order ^/P, thermalization can be confirmed. That means that 
macroscopic fluctuations of the solutions decay at the benefit of microscopic fluctuations. 

The mathematical analysis for the microkinetically regularized wave equation consists in 
two parts. First we present some analytical and numerical results on the propagation of 
phase boundaries and thermalization effects. Despite the fact that model is conservative, 
it exhibits the hysteretic behavior. Such a behavior is usually interpreted in macroscopic 
models in terms of dissipative threshold which the driving force has to overcome to ensure 
that the phase transformation proceeds. The threshold value depends on the amount of the 
transformed phase as it is observed in known experiments. 

Secondly we investigate the dynamics of oscillatory solutions. Our mathematical tool are 
Young measures, which describe the one-point statistics of the fluctuations. We present a 
formalism which allows us to describe the effective dynamics of rapidly fluctuating solutions. 
The extended system has nontrivial equilibria which are only visible when oscillatory solutions 
are considered. The new method enables us to derive a numerical scheme for oscillatory 
solutions based on particle methods. 




Figure 1 . Nonmonotone stress strain relation and nonconvex elastic energy 

1. Introduction 

In this work we study the onset and propagation of oscillations in a model for martensitic 
phase transformations. Consider the simple tension of a rod of nonlinear elastic material 
with a diagram a = shown in Figure 1, where a and e are tensile stress and strain, 
respectively. The diagram shown in Figure |l| corresponds to a nonconvex free energy W{e) 
and is typical for materials with martensitic phase transformations (PT). We will label 
the phases corresponding to the branches BC and EF (-) and (+) respectively; branch 
CDE exhibits an unstable intermediate state. The principle of minimum of free energy for 
such materials results in formation of a phase mixture and in a macroscopic deformation 
along the line BDF at a constant stress ctm (Maxwell line), rather than along the lines 
BCDEF. Considering a quasi-static experiment at some stress ctq, a macroscopic portion 
of the material is deformed from e~ to . At the next total strain increment Ae, the 
next portion of material undergoes PT. During the whole process of PT, the strains in 
the phases (-) and (+) are fixed and equal to e~ and the stress is constant a = 
and Ae is related to the increment Ac of the volume fraction of a new phase. 

It is clear, that at quasi-static treatments we do not follow the original stress-strain 
curve e <y{e), but we rather introduce a strain jump from e~ to at constant stress 
for each transforming particle. To be consistent with the original stress-strain relation, 
the concept of fluctuating stresses was introduced in 0. Fluctuating stress Ufluct is 
defined as cTfluct = cri^e) — ctm in the interval [£",£"*"], it varies in accordance with the line 
BCDEF, relative to the line BDF. This fluctuating stress overcomes the energy barrier 
BCD, but due to the very small time interval of its action, it is not taken into account 
in the continuum relations. In the equilibrium process, the time averaged value of cf^^ict 
is equal to zero. 

This concept is useful for the development of some kinetic models |]T], ^ which only use 
some important properties of the stress-strain curve a = a{e) but not the precise course 
of the unstable part. 

One of the main unsolved problems is to find a proper physical model for fluctuations 
and to include it in the continuum description of PT. A possible way is to take into 
account microinertia in constitutive equations, i.e. to consider a constitutive equation 



(1.1) 



am = cr{e) + (3e, 
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where cxm is the macroscopic stress. The constitutive assumption ( |1 . 1| ) and the conserva- 
tion of hnear momentum lead us to the following evolution equation: 



(1.2) iM = d,{a{d,u) + pd^il), X G (0, 1), p > 0, /? > 0, 

where u : x ^ M is the deformation of an elastic body i7 = (0, 1) C M at time t and p 
is the mass density. The coefficient (3 describes the inertia of microscopic fluctuations, for a 
derivation see Chapter ^. For reasons which will become clear there, we will refer to system 
( |1.2| ) as microkinetically regularized wave equation. If appropriate boundary conditions 
are imposed, the generalized energy Jq{^u'^ + W{dxu) + ^{dxu)'^) dx is conserved. The 
aim of this paper is to present the first mathematical results for such a model. 

Equation ( |1.2|) constitutes a well posed Cauchy problem. In the following section we 
present a continuum mechanical derivation of the microkinetically regularized wave equa- 
tion based on an approach developed in ^, Q to describe the dynamics of phase trans- 
formations. Since in the derivation finite wavelengths are assumed, it is not possible to 
show that ( |1.2| ) approximates the phase transformation problem if the strain oscillates 
at a wavelength of order However the obviously very interesting behavior of the 

solutions motivates us to analyze the dynamics of oscillations with modern mathematical 
tools like Young measures. Many of our results are also valid for the improved Boussinesq 
equation 

U Uxx ~l~ )xx ~^ ^xxi 

see Chapter |^ for references. 

In Chapter ^ important phenomena of the dynamics are investigated numerically and 
analytically. We show that the system possesses solitary-wave solutions. In the numerical 
simulations quasi-dissipative behavior can be observed. Due to the nonmonotonicity of 
a the solutions exhibit high frequency oscillations in space after times, which are pro- 
portional to a/5 even if the initial state is smooth. These spatial fiuctuations can be 
interpreted as heat. The introductory problem motivates the development of mathemati- 
cal tools which simplify the study of oscillatory solutions. We expect that the results will 
be useful for future work. We find a critical threshold value for the minimal driving force 
of the transformation process. The threshold value depends on the amount of transformed 
energy. 

There exists a rich literature concerning regularizations of quasilinear wave-equations 
with nonmonotone stress-strain relation, see e.g. for a starting point. The problem 
lies in the ill-posedness of the Cauchy problem for the quasilinear wave-equation pu = 
dx(y{dxU). The hope is that the regularization selects the physically correct solution if the 
strength of the regularization tends to 0. Another approach consists in the assumption 
of a kinetic relation to determine the velocity of a moving phase boundary. Using this 
additional constitutive law one can show that the Riemann problem has a unique solution, 
see [0, H . Indeed the solutions of the Riemann problem with viscosity and Van-der-Waals 
forces generate the kinetic relation, 

We are not so much interested in the singular limit /? — but more in the qualitative 
behavior of the solutions if (3 is nonzero. The system distinguishes in a natural way 
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between two different scales (microscale and macroscale). This phenomenon allows us to 
characterize the dynamics of microstructure in a unique fashion. To this end we introduce 
the notion of Young measures which can be viewed as generalized limits of sequences of 
highly oscillatory functions. One major disadvantage of Young measures is that they do 
not work well together with pseudo differential operators of degree like the Helmholtz 
projection which maps a vector field onto the gradient part. But this is essential if space 
dimensions greater than one are considered. Therefore the analysis is restricted to the 
one dimensional case. 

Young measures have been used successfully for phase transformation problems. In [ p!0|] 
they are applied to the three dimensional crystallographic theory. Using methods from 
calculus of variations Kinderlehrer and Pedregal showed in [|lT] that dynamical 
models for phase transformations admit Young measure solutions. The approach outlined 
in this work has been applied before to the viscoelastic system with nonmonotone stress- 
strain relation pu = dx{<j{dxu) + Pd^u), see [|T^ . 

The mathematical analysis for oscillatory solutions of the microkinetic model consists 
of two parts: short-time behavior and long-time behavior. 

In Chapter ^ we derive an effective equation for Young-measure solutions. These gen- 
eralized solutions are unique and they approximate the dynamics of highly oscillatory 
solutions over finite times. This approach can be viewed as an analogy to homogenization 
theory. There the aim is the derivation of a macroscopic law if the the material properties 
are oscillatory. In our case it turns out that the dynamics of oscillatory solutions can be 
reduced to a transport equation which has structural similarity to the Vlasov-Poisson sys- 
tem. The Vlasov-Poisson system describes the evolution of electrically charged plasmas 
(see e.g. 0). 

The microkinetic system exhibits some kind of ergodic behavior. Even if we start with 
a smooth initial condition the solutions soon become oscillatory, cf. Figure |^. Especially 
in the region x G [0.1, 0.2] the solution looks very much like a Young measure. 

In Chapter 5 we demonstrate that the extended system has a large number or equilibria 
which are not generalized limits of sequences of stationary solutions of ( |1.2| ). The physical 
interpretation of this result is as follows. Despite the very dynamical oscillatory behav- 
ior of each material point, macroscopic statistical characteristics, like volume fraction of 
new phase are stationary for equilibrium solutions of the extended system. Thus Young 
measure are an important concept for the finding of macroscopic information using non- 
stationary oscillatory solutions. This observation has been made earlier for the equations 
of ideal incompressible magneto hydrodynamics, see |T^. There the authors analyze 
so called relaxed equilibrium states which correspond to highly oscillatory nonstationary 
solutions of the evolution equations. 

If we consider the nonphysical limit p — it can be shown, that almost every solution 
tends to a stationary Young measure as time goes to infinity. This is a genuinely infinite 
dimensional result since the system is Hamiltonian. In finite dimensional Hamiltonian 
systems the phase-space volume is conserved, thus generic solutions do not converge to 
an equilibrium. 
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In the long-time behavior of oscillatory solutions has been investigated for the 
viscoelastically damped system. There it is found that the microstructure region shrinks 
when t tends to infinity. This contrasts to the microkinetic system where microstructure 
can be generated in the long-time limit. 

The convergence towards an equilibrium Young measure can be understood as thermal- 
ization, i.e. the fluctuations which are at the beginning only visible at the macroscopic 
level are transformed into microscopic fluctuations. Hence on the macroscopic level the 
solution becomes stationary. This behavior does not correspond to the numerical results 
of E. Fermi, J. Pasta and S. Ulam in |]T6[ where a discretized nonlinear wave equation 

pEi = a{ei+i) - (r{ei-i) 

is investigated. They observed a recurrence phenomenon, i.e. after fixed time intervals 
the solution returns very close to the initial deformation. 

The study of the microkinetic system is concluded by a systematic numerical investi- 
gation of the thermalization behavior in the case p > 0. In Chapter]^ a numerical scheme 
based on particle methods is developed and simulation results concerning the long-time 
behavior of solutions for values of p greater than are presented. It turns out that the 
spatial coupling due to the nonvanishing of p complexifies the long-time behavior con- 
siderably. The solutions develop microstructure, but there are no indications that the 
solutions tend toward a stationary Young measure if p > 0. However, a major portion of 
the initial macroscopic kinetic energy is transformed into microscopic fluctuations. 

2. Constitutive equation for material with microscopic inertia 

For a formal derivation of the constitutive equation we consider the c?-dimensional 
setting, d = 1,2,3, although the subsequent analysis is restricted to the one dimensional 
case. For simplicity geometrically linear elasticity theory is used, we insist however that 
the derivation is completely analogous if geometrically nonlinear theory is used. In this 
case we can use the deformation gradient and the first nonsymmetric Piola-Kirchhoff 
stress tensor as conjugate variables. For the convenience of readers which are not familiar 
with the notations in continuum mechanics we repeat the basic definitions. Let f2 C M'^ be 
the reference configuration of an elastically deformable body. The superscript t denotes 
transposition and V is the gradient operator. We use bold letters to indicate tensors and 
vectors, the coefficients are written with normal letters. The dot " ■ " and " : " denotes 
the contraction of tensors over one and two indices: (A • B)ik = ^ijBjk = AijBjk and 
A : B = tr(A ■ B^) = AijBij, where Einstein's sum convention is adopted for the last 
identities. By A® B we denote the fourth order tensor with components (AijBki). 

Let us consider the representative volume u) G Q of the microheterogeneous material 

bounded by a surface du with the unit normal vector n . We assume that u is much 

smaller than Q and that the center of mass is located at the origin. Let r be the 

particle position, u the continuous displacement field, i.e. for every material point r in 

the reference configuration the position at time t is given by u{t,r). By a- we denote 

the stress tensor and £ = Vu is the deformation gradient tensor, where ~ means the 
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local value of parameters as opposed to the macroscopic counterpart. Let us introduce in 
the usual way the macroscopic stress cr^ and deformation gradient s.^ tensors 



(2.1) cr^ = vol(u;) W r (g) cr • n dS ; = vol(cj)~W u(g)ndS. 

J did J duj 

The macroscopic stress tensor is the critical quantity which allows to analyze the macro- 



scopic behavior. Definition ( |2.1| ) can be motivated as follows. Let be the total energy 
contained in the representative volume, then the power of the external forces can be 
computed with the formula 

(2.2) ^E^= f (5r*-S)-nd^. 

Assuming that the velocity distribution corresponds to a macroscopically homogeneous 
deformation on the boundary, we find that n(r) = uq + Sm ■ r holds for all r G du. 
Plugging the Ansatz in (p.2|) one finds that 



/ (cr* ■ {srn ■r))-ndS+ {cr ■ n) ■ iio dS 



=e\^ : / r (g) cr ■ n dS* + tio • / (cr ■ n) dS* 

J duj J duj 

where Fm = /g^(^ ■ ''^) dS* is the macroscopic force. Thus in absence of macroscopic 
translations the power of the external stresses, through which the representative volume 
interacts with the neighborhood, is equal to cr^ : £m, as for material points. The symme- 
try of the stress tensor in the geometrically linear theory is corresponds to the symmetry 
of the Cauchy stress tensor in the geometrically nonlinear theory. With respect to both 
theories the macroscopic stress tensor is in general not symmetric, it is possible to transfer 
energy by rotations. The macroscopic forces will be neglected in the further considera- 
tions. Using the Gauss theorem and the equations of motion without body forces 

(2.3) V-d- = pit , 

where p is the mass density, we obtain from eqs. (|2.1|) and ( p.3|) 

(2.4) cr^ = (or) + {pit (g) r) ; = (e) , 

where (...) = vo\{uj)^^ J^{...) dx is the averaging over the volume u. When the macro- 
scopic stress tensor is approximated by the averaged stress tensor the resulting errors can 
be expressed as accelerations or inertial stresses. The averaged value (ct) neglects rapid 
oscillations in the interior of the representative volume. The kinetic energy of these oscil- 
lations is responsible for the inertial stresses. It is crucial that the inertial stresses have 
a regularizing effect. They enlarge the inertia of the material in such a way that a well 
posed Cauchy problem can be derived even if for the averaged stress tensor a nonmonotone 
dependency of the macroscopic deformation gradient is assumed. 

If the diameter of u tends to then in (B~3!) the second term vanishes and 



(2.5) CTm = {^) 
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i.e. the macroscopic stress cr^ is equal to microscopic stresses averaged over the infini- 
tesimal volume cr^ := {a-) . The same holds if the body is in an equilibrium state. The 
next step consists in prescribing the constitutive equations for a (in the micromechani- 
cal approach) or directly for cr^ at phenomenological consideration. In order to describe 
martensitic PT in elastic materials a nonconvex elastic potential is used, which results 
in nonmonotone stress-strain relation cr^ = f {Sm) ■ Due to the ill posedness of the 
dynamic boundary-value problem for such type of materials, a number of mathematical 
and physical regularization methods have been suggested. One of them consists in the 
introduction of viscosity, i.e. in considering a viscoelastic material. In the given paper 
we would like to limit ourselves to nondissipative materials. Another approach is related 
to introduction of characteristic material size, e.g. by introduction of higher strain gra- 
dients V , V V £m , . . . or nonlocal stress-strain relations. As it is shown in [|I| , in 
a framework of micromechanical consideration, the macroscopic deformation gradient B 
and the work-conjugated third-order tensors of the hyperstress M can be defined by the 
formulae 

(2.6) B = Yo\{uj)-^ f e®ndS] M = ^"^^"^^ / r0r®a-ndS. 

With the help of the Gauss theorem and the equations of motion for static case we obtain 

(2.7) B = (VVw) = (Vi) ; M={r®a). 

The tensor M is nonzero only if the diameter of the representative volume u is finite, this 
follows from equation ( |2.7] )2 by considering the limit diam(co') 0. At macroscopically 
homogeneous strains u = uq + Sm' r on duj we have e = Sm on the surface du), and 
according to equation (|2.6| ) B = . li the macroscopic stresses on du are homogeneous 
a-n = (Tm-n holds and according to equation ( PT7| ) M = (using (r) = 0). Conse- 
quently, the tensors B and M characterize the heterogeneous distribution of strains and 
stresses inside the representative volume. Both vanish for macroscopically homogeneous 
conditions even if the representative volume is finite. 

At the same time for finite representative volume, the inertial stresses cr, := (pr^u) 
are nonzero for macroscopically homogeneous stress state, but high accelerations it . As 
PT represent a highly dynamic process, the contribution of inertial stresses can be in fact 
very important for finite representative volume. 

Our next aim is to approximate the compact integral operator {pr®u) by an unbounded 
differential operator. Expanding u around to the first order we have u{r) = u{0) + 
£(0) ■ r + 0{diam{ujf\\\/^iL\\L--(n))- Thus 

{pr ®u{t,-))ij = ^^^l / XiUj{0)dx + '^ / ejk{t,0)xiXkdx> + 0{dieim{Ljf\\V'^u\\Lo 

. k=i ^'^ J 



=0 



d » 

--P^£jk{t,0)^^ / XjXfcdx + C(diam(a;)^||V^w||L->(f7))- 
k=i -''^ 
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Due to the hnearity of the integral operator the estimates are stable under weak conver- 
gence, i.e. we may replace e by Dropping the error term we obtain 

(2.8) cr^ = cr„ + (pr (g) r) • 

It is clear that the approximation looses validity if the oscillations of 'im occur on a length 
scale of order diam(cij). 

The final step in our derivation consists in taking equation ( p.4| ) as a definition for 
CTm- This introduces a correction term, which ensures that the averaged stress contains 
enough information to predict the macroscopic dynamics. This new stress is conservative 
and related to the acceleration, therefore we refer to it as microinertial stress. It has a 
similar effect as a viscoelastic damping, namely the spatial regularization of the velocity. 

At one dimensional tension, when only tensile stress a and strain e are important, 
equation (|2.8|) results in 



(2.9) a^ = a, + Pe- (3 = {pr') , 

where the coefficient (3 characterizes the inertia of the representative volume at the fixed 
center of mass, i.e. microscopic inertia. Equation ( |2.9| ) resembles to a Lagrange equation 
in mechanics of discrete systems with the generalized coordinate e under the action 
of external and conservative generalized forces. For this reason we will denote 
the kinetic energy of the representative volume | /5 as microkinetic energy. If the 
representative volume is a cube with side length 2a and the mass density is homogeneous, 
then /? = f a2. 

Equation ( |1.2D has been derived in the case a{e) = Ee {E is the Young modulus) to 
describe longitudinal waves in a slender elastic rod [0, p. 428], the equation is known 
as Love's modified wave equation. The inertia of radial deformations are taken into 
account; this allows in contrast to the classical one- dimensional wave equation to describe 
wave propagation problems in rods where the size of the diameter is comparable to the 
wavelength. In the case a{e) = e + e^, p = 2 the equation is known as improved Boussi- 



NESQ equation, it has been derived in [|T8| to describe ion-sound waves. The classical 
Boussinesq equation 



(2.10) il = d;{u + + d;u) 

which has been derived in 1872 to describe shallow water waves [|l^ has the shortcoming 
that the Cauchy problem is ill posed. Thus it cannot be used to analyze wave propagation 
problems numerically. The dispersion relation of both Boussinesq equations 

tu^ = /c^(l — k') (classical), 
cu^ = A;^(l + k')~^ (improved). 



are equivalent to the fourth order for \k\ -C 1, this motivates the nomenclature. We will 
focus on the short wave-length limit \k\ oo. In this limit the improved and the classical 
Boussinesq equation are not related. 
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For odd numbers p the mapping e — *■ a{e) is monotone, for even numbers the solutions 
can blowup in finite times, see, e.g. [^. We are interested in the case that a is non- 
monotone but the stored elastic energy function W (W' = a) is coercive, i.e. W{e) oo 
as £ — *• oo. To our best knowledge this case has not been studied in the literature yet. 

Note, that an inertial term in stress-strain relations was introduced and made consis- 



tent with thermodynamics in [Q, 3. Earlier Valanis [^] suggested inertial contribution 
in constitutive equations for internal variables, and a way to satisfy the second law of 
thermodynamics in this case. A general theory, which includes generalized inertial forces 
for each thermodynamic variable, was developed in 0, |^ . It was shown, that in particular 
cases thermodynamic inertial forces were well known for thermal conductivity, moisture 
propagation in a colloid capillary body and turbulent transfer, but the equations were not 
consistent with the second law of thermodynamics. 

3. First discussion of the equations of motion, solitary waves 

In this section we state the precise evolution equations and explain central proper- 
ties. To motivate the consideration of oscillatory solutions, simulation results of special 
solutions are presented. 

We set i7 = (0, 1) C M. Using the constitutive equation (|2.9|) where the subscript of cXm 
is dropped, we obtain the following nonlinear partial differential equation which governs 
the time evolution of u: 

(3.1) pu = d^{cr{d,u) + f3d,u), /3 > 0, p > 0, 

u\t=Q = Mo, u\t=o = vo (initial conditions), 

u\x=o = 0, {a{dxu) + pdxu)\x=i = (Neumann-boundary condition). 

The chosen boundary conditions model an experiment where the rod is free at the right 
boundary and fixed at the left boundary in order to exclude translations. 
We assume that a{-) is a nonmonotone function with three zeros e~ < < which 
satisfies lim|e|^oo£^ ■ (^{e) = 00. For mathematical simplicity we do not consider boundary 
forces to control the experiment. Instead the shape of the mapping e 1—^ a{e) is varied 
with a scalar parameter /i in order to stabilize or to destabilize a certain phase, therefore 
no assumption that J^__ a{e) de: = is made. The nonconvex function W{e) = (7(7) d7 
is usually called double-well potential, the two local minima correspond to two differ- 
ent phases which are stable with respect to small perturbations. Global existence and 
uniqueness of solutions to (|3.1| ) is clear, although a formal proof is not provided in this 
work. 

If /5 = system ( p.l| ) describes the purely elastic behavior of the rod where no microki- 
netic forces are taken into account. It is clear that the initial-boundary value problem 
is well posed only if a is monotone and thus equation (|3.1|) is hyperbolic. By assuming 
a nonmonotone stress-strain relation we have introduced a Hadamard-instability which 
causes the solutions to start oscillating heavily and blow up in finite time. The term 
Pdlu is of higher order thus for /? > global existence and uniqueness of solutions can be 
established. 
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System ( p.l|) is the evolution equation associated to a Hamiltonian system. The Hamil- 
tonian (energy) which is conserved by the dynamics is given by 

H{q, p) = I 7:P ■ Bp + W{q) dx, 

where g(t,x) = d^uit.x) E L'^{n), p{t,x) = {p - (3dl)u{t,x) G H-\n) and Bp E H^Q) 
is the solution of the elliptic equation 

{p-(3dl)v = p, (9^ti|3.=o = 0, f 1^=1 = 

e.g. V = Bp. The symplectic structure is given by the skew-symmetric operator J = 
. In u, M-coordinates the energy takes the form 









E{u, u)= I (f + W{d^u) + ^{d^iif) dx. 



The expression -f^mac = f /^g^ dx is called macrosco])ic kinetic energy, E = j^^^ W{dxu) dx 
elastic energy and -fTmic = f j'x&n^^^'^)'^ '^■"^ 'microscopic kinetic energy. 

3.1. Qualitative and quantitative properties of typical solutions. The term I3d1u 
in ( p.l| ) models the inertia of microscopic structures, e.g. phase boundaries. To give an 
illustration, we consider a phase boundary which is at time i = at the position x = 
and moves to the right with velocity c. This corresponds to the initial value 

1^ X — k; X E [ —2c; x E 

The variable n <^ 1 describes the thickness of the transition layer between the two phases. 
With respect to the stored energy function W{e) = ^^e"^ — 1)^ we obtain the energy 

E{uo, vo) = (2 - Ik)pc' + + 2^. 

This example shows heuristically that the velocity c of the phase boundary has vanishing 
influence on the macroscopic kinetic energy and the elastic energy, if k is close to 0. Only 
the microscopic kinetic energy is able to measure the velocity of the phase boundary. 

The solution of system ( |3.1| ) generated by the above initial condition is by no means a 
traveling front in the variable dxu{t, ■) as it can be seen in Figure ^ where the graph of 
dxu(t, ■) is plotted for different values of t. We have used a{e) = E{e'^ — e), the Young 
modulus E = 1 MPa and the density p = in the numerical simulation, for a precise 
description of our numerical scheme, error estimates and additional results on the long- 
time behavior see Chapter ^. The solution depicted in Figure ^ has been computed with 
mirrored boundary conditions i.e. Dirichlet-boundary condition at x = 1 and Neumann 
at X = 0. It is obvious that the transition region with thickness k generates a strongly 
oscillatory region which is located in the interval [0, 0.2]. In addition one can identify a 
pulse at X = 0.45 which is moving at constant velocity to the right. 

Solitary waves have been found and analyzed for modifications of the Boussinesq equa- 
tion ( [^.lOD as well as for the improved Boussinesq equation. Formal calculations show 
that the former equation is completely integrable and the latter is not since it does not 
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Figure 2. Initial value 9^mo(-) 



dxU at time t 



0.4 (p = 1, /5 = 10" 




2, K 



0.1 and a snapshot of the strain 



pass the Painleve test. This implies that solitary waves may exist but the interaction 
of two pulses is not elastic, i.e. energy is transferred into uncorrelated fluctuations, see 
The analysis in p2| shows that also for the microkinetically regularized wave 
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equation it cannot be expected that it is completely integrable, this is confirmed by the 
results of the numerical simulations. 

We find two kinds of solitary waves which we denote by nudeation pulse and elastic 
pulse. The amplitudes of the nucleation pulses are bounded from below, they traverse 
the nonmonotone sector of a completely. The amplitude of the elastic pulses may be 
arbitrarily small. Elastic pulses have nothing to do with the nonmonotonicity of a, they 
are closely related to the solitary waves which have been analyzed in the Boussinesq 
equation. 

3.2. Nucleation waves. To discuss traveling waves we make the usual Ansatz (9i.M()f:, x) = 
p{x — ct) where p{-) is the initial strain and c is the velocity of the propagation. Plugging 
p into ( |3.1| ) we obtain the ordinary differential equation 

pcV = a{p)" + pcY". 

Integrating the above equation twice yields 

(3.3) pc^p = a{p) + /?cV + A + Ax, 

where A and A are the constants of integration. We are interested in solutions p satisfying 
lim|a,|_+oop(x) = e^, i.e. nucleations within the phase which is given by the local minimum 



e = oi W, cf. Figure |I[ Using equation 
we end up with the second order equation 



we find that A = pc^e^ and A = and 



(3.4) 

The solutions of 



Q 2 II ( \ \ 2 2 + 

pc p = —cr{p) + pc p — pc e . 
move along the equipotential lines 
\Pc^{j)'f + W{p) - \pc^p^ + pc^e+p = const. 



In order to find a homoclinic orbit attached to the point y'^ = {e'^, 0) G the point 
has to be a local maximum of the function {e i— > W{e)—pc'^e{^e—e^)), hence pc^ > <j'{e^) 
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0.2 04 06 0.8 1 

Figure 3. The shape of elastic and nucleation waves for p = 1, /i = 0, 
c = 2 and /5 = 0.025 (-), 0.0125 (■ ■ ■ ) and 0.0025 (-) 

The amphtude A of the pulse is |£:+ — 7I where 7 is the smallest respectively greatest zero 
of 

(3.5) ^{e'Y + W{e) - ^pc'e^ + pch+e - W{e+) - ^pc\e+f. 

For the choice W{e) = — 1)^ + p{£ — which has been used in the numerical 
simulations, the shape of the pulse can be computed explicitly: 

p{x) = 1 + H \ 



1 - f ± ^ipc2 + lp{p + 3) cosh 



pc2-2(l-M) 

I3c- 



,2 X 



The role of p will be explained in the following section. Thus we have obtained two 
families of pulse solutions (pei and Pnuci) moving at velocity c. If p = one finds the 
simple velocity-amplitude relation A{c) = max^^^ — p{x)\ = Cy/2p ± 2. 

For p = 1, c = -\/2, p = the amplitude A = 4 is approximated well by the numerically 
computed pulse. Note, that in the case p = pulse solutions satisfying lim^^^^^p^x) = 
do not exist. Of course these results also hold true with minor modifications if the stored 
elastic energy W is given by a more general function than a polynomial of fourth order. 

3.3. Propagation of a phase transformation front. To discuss the simulation results 
concerning the considerations in the introduction, we would like to study the propagation 
of a phase transformation front. From the preceding discussion of the pulse solutions 
it is clear that front solution can not exist. Front-like solutions require that the energy 
W{e) — ^pc^e^ + \e (cf. ( p.3|) ) has two local maxima of same height at the positions e = e~ 
and e = , but this is only possible in degenerate cases which we will not consider in the 
present work. 

The model is intended to describe the dynamics of phase transformations when no dis- 
sipation is present. Since the solutions do not have the possibility to dissipate kinetic 
energy, the fluctuations due to the phase transformation remain present after the trans- 
formation region has propagated. Therefore we expect to observe wild oscillations in the 
wake of the nucleation pulse. These oscillations should be interpreted as heat since they 
do not give a contribution to the kinetic energy but only to the microkinetic energy. 

We consider a scenario where phase 2 is slightly destabilized i.e. the boundary force 

S is below the Maxwell-line. It is unpleasant that the local minima of the Gibbs' energy 
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W{e) — S ■ £ are only given implicitly, this makes the algebraic manipulations intricate. 
Therefore we model the effect of S 7^ by taking a{6) = —e + fi- (1 — e^) where yU G M is 
a control parameter, which allows us to destabilize (/i > 0) or to stabilize (yU < 0) phase 2. 
For <C 1 we obtain the relation 

3. 

~ — 



The simulation of the solution defined by the initial datum ( |3.2D shows exactly the ex- 
pected behavior, see Figure ^. One can observe a phase transformation front characterized 




X 




X 



0.2 0.4 0.6 0.8 1 
Figure 4. Two snapshots at time t = 2 and t 
phase-transformation region gets stuck at Xte 
/i = 0.5). 



0.2 0.4 0.6 0.8 1 

= 3 show, how a propagating 
„ = 0.6 (p = 1, /? = 10-2, 



by the the tip of the pulse at x = 0.45 moving from left to right through the material, 
leaving a wake of oscillations behind. It is interesting that the phase transformation 
front comes to a stop when the gained energy is transformed too quickly into microscopic 
fluctuations. This obstructs the activation of the stable phase and the phase transfor- 
mation may terminate, see Figure ^ This astounding effect demonstrates that the study 
of the qualitative dynamic behavior requires a good understanding of the high frequency 
oscillations in space. If the sign of fi is reversed, phase 2 is stabilized and the phase 
transformation does not even start, i.e. hysteresis is observed. It is very remarkable that 
allowance for microkinetic energy results in this phenomenon. In macroscopic models it 
is interpreted in terms of dissipative threshold which the driving force has to overcome 



in order to phase transformation proceeds (see e.g. ||23| , |2J]). Moreover in the numerical 
simulations this threshold depends on volume fraction of transformed phase, which cor- 
responds to macroscopic description of known experimental results [p5| , In Figure |] 
the position where the phase transformations stopped is plotted over different values of 



4. The propagation of oscillations 

The Hadamard instability which is present in equation ( |1.2| ) if /5 = is not completely 
removed by the regularization term pdlu. The onset of the fluctuations is delayed, we 
expect that after time T the oscillations will occur on a length scale h = ^f^jT . Their 
presence makes a rigorous discussion of the qualitative dynamical behavior very difficult. 
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Figure 5. Position of the stuck PT- front for different values of fi. 



Our approach is to study the oscillations which are brought via oscillatory initial states 
into the system. We do not let P tend to 0, instead we consider for fixed /5 > a sequence 
of oscillatory initial states and ask: "How do the oscillations propagate as the time in- 
creases?". It turns out that they remain stationary in Q, but the stochastic properties 
change in a well defined fashion. To give a motivation one could think of taking for /3 ^ 1 
the solution at time t = 1 (which should look very similar to the simulation result depicted 
in Figure 0) as a new initial value. 

In the derivation of the constitutive law (|2.9|) we have assumed that local stresses and 
strains a and e can be replaced by the homogenized counterparts am and Em- Thus 
the constitutive law ( p.9[ ) looses its validity if oscillations occur which have a wavelength 
shorter than ^/J3, the size of a characteristic volume. On the other hand there are only 
few systems known, which support oscillatory solutions. The authors are not aware of 
any model where the dynamic creation of microstructure can be studied rigorously. We 
will ignore the modeling restriction and concentrate on the discussion of of oscillatory 
solutions to the partial differential equation ( ^TII) . 

The mathematical analysis of oscillatory solutions consists in several steps. First an 
integro-differential equation, which allows a clearer view to the structure of the equation 
is derived. Then we define Young measures and oscillatory solutions in a mathematical 
rigorous way. The propagation of oscillations is characterized by the extended system 
( [4.2|) , see Theorem 4^. In Chapter |^ the stationary solutions of the extended system are 
characterized and discussed. If p = the long-time behavior or the microkinetic system 
can be understood completely. Almost every solution generates a nontrivial stationary 
Young measure as time tends to infinity. The deformation u{t) becomes stationary but 
the long-time limit is not an equilibrium solution of (|1.2| ) (Theorem ^.3|) . 

In order to simplify the analysis for system ( p.l| ) we transform it into an integro- 
differential equation which does not contain spatial derivatives. The new variables are 
given by {e,ri) = T{u,u) = {dxU,dxu)- By differentiating equation ( ^.1| ) once by x one 
obtains 



{p-(3dl)e = dla{e). 
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Using the operator B which has been defined in section 1^ we find that 



e - e\r,=i = pB ( dl-cr{e) - e\^=i 



P 

(p-/39^)i(a(5)-a(5U)) 



=0 



Thus {e, 7]) satisfy the following integro-differential equation 

(4.1a) £ = r/, e\t=o = Sq e L^{Q), 

(4.1b) T] = -|(a(£) - k), r]\t=o = Voe L^iil), 

(4.1c) k = pBa{e), 

which defines us a dynamics in L^(f2,R^). The evolution preserves the regularity of the 
solutions i.e. if {eo,Vo) e C\Q,R^) then {e,ri){t) G Ci(fi,M2) for all t > 0. The most 
important property of this system is that the equations for e and rj are hyperbolic with 
wave speed 0, thus the solutions cannot develop shocks. This transport equation allows 
the system to carry solutions oscillating rapidly in e and rj. Since the spatial coupling via 
pB is compact, k cannot fiuctuate in space. 

4.1. Definition of Young measures and oscillatory solutions. A relevant tool for 
the study of oscillations are Young measures (or parameterized measures) which have 
gained of importance during the last twenty years. A Young measure \E' G Y(f2,]R'^), 
(i G N, is a mapping from Q into the set of probability measures on M"', which is denoted 
by PM[M.'^). For each x & Q the probability density '^{x){-) measures the probability 
that the strain and the strain velocity {dxU^"'\ dxU^"'^) take the value {s,r]). 

In the context of oscillations Young measures are conceived as a family of one point 
distributions which are usually generated by a sequence of rapidly oscillating functions. 
In order to clarify the term "generated" we have to become a little bit more technical and 
introduce notions from probability theory. Let /'•"^ : — > M be a bounded sequence of 
functions. For a given point a;o G ^2 we choose a positive number r, a positive integer n 
and a step length h. We intend to consider the limits h —>■ 0, n —>■ (yo and r ^ in this 
order. The triple (r, n, h) defines a finite sample consisting of the values f^"'\xo + hk)^ k G 
{q G Z|a;o + 'qh & [xq — r, + r]}. Now we let h tend to 0. This generates a unique family 
of probability measures i/(a;o, r, ra, ■) G PM{W^) which is parameterized over (xojT, n) G 
VL X X N. If for almost every xq the limit z/(a;o) = limr^o linin^oo ^{xq, r, n) exists 
(in the vague sense) then we say that /^"^ generates the Young measure u : VL ^ PM(M'^) : 
Xo ^^ z/(xo). 

An alternative characterization, which is more suitable in the context of partial differ- 
ential equations, is given by the duality of probability measures and continuous functions. 
A sequence f^'^^ generates a Young measure \E' if and only if for every continuous function 
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g : R'^ 



lim - {^{x),g)^^^, weakly in L\Q). 



By "(-, ■)" the duality between probability measures and continuous functions is denoted. 
To obtain a clear distinction between the variables the integration variable is used as an 
index, i.e. 



{^{^),9)yeR^ = I g{y)-'^{x,y)dy 



if the probability measure \l'(x) has a density y). By slight abuse of notation we will 
adhere to this function notation even if the measure is singular. It is our aim to emphasize 
the geometric viewpoint by this notation. 

We say that the sequence of solutions m*^"-* defines an oscillatory solution if the initial 
values {dxU'Q\dxV^^) generate a nontrivial Young measure \E'o- 

4.2. Effective evolution equations. Spatial oscillations (microstructure) are not damped 
away by the dynamics of the microkinetic model. A possible characterization of the time 
evolution of oscillations is given by the following problem: 



Given a sequence of initial states to the microkinetic system (ug"'', fg"'') so that 



X] 



{dxUQ , O^Vq ) generates a unique Young measure \i/o as n — oo. Does 
{dxU^"'\dxU^'^^){t,x) generate also for t > a unique time-dependent Young measure 



It can be shown that the answer is positive, the time depending Young measure "^{t, x, y) 
is a distributional solution of the transport equation 

(4.2a) dt^{t, X, y) + Vy {^{t, x, y) ■ (2/2, -^(^(l/i) - k{t. ^)))*} = 0, ^|t=o = ^0, 
(4.2b) (p - Pdl)k{t, x)=p {^(t, X, y), a{y^))y^^, , 9,fcU=o = 0, kl=i = 0. 



The transport term is exactly the right-hand side of system ( [4.1|) . Equation ( [4.2a| ) has to 
be read in the weak sense, i.e. for every test function g G C^(]R^) 

(4.3) ^ (vj/, g) - (vl/(t, X, y), V g{y) ■ (2/2, -^(^(z/i) - ^(t, x)))*) = 



for every x G fi. System (|4.2| ) is a extension of (p.l|), i.e. every solution of (|3.1| ) is a 
special solution of (|4.3| ). We formulate this fundamental assertion in a 



Proposition 4.1. Let u{t,x) be a solution of (|3.1| ) and e = d^u, rj = d^u- Then the 
Young measure 

^'(t,x) = 5(^e,r)){t,x) 

where 6y is the Dirac distribution with mass aty eM."^, is a distributional solution o/( |4.2| ). 
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Proof. Since the equivalence of and (|4.1| ) has already been established, we set 

y{t,x) = T{u{t,x),u{t,x)). By inserting Sy(^t,x) into ( ^l3|) we obtain 

^^9{y{t)) - Vgivit)) . {y,, -i((7(yi) - k{t, x)))' = 

by (O)- □ 



It can be shown under reasonable assumptions that the solutions of ( ^4.2| ) are unique and 
exist globally in time. The central assertion which justifies ( [4 .21) as an effective equation 
for oscillatory solutions is the following continuity result. If the oscillations of the initial 
state can be captured with a Young measure it suffices to take an initial value for 
( [4.2|) . The unique solution at time t is a Young measure which describes the oscillations 
at time t. 

Theorem 4.2. Let {dxU^f^\ d^v^Q^) G L^(fi, M?) generate the Young measure G Y(n, R^) 
and u^^\t) be a sequence of solutions to the microkinetically regularized wave equation 
so that ((9^M("),9^M("))(t = 0) = Then generates 

a Young measure ^E'(t), which is the unique solution of ( [4.2|) satisfying the initial condition 
^(0) = ^0- 



Proof. For a proof see |26|, ^7j. □ 



It is not at all clear that Young measures are always the right objects to capture oscillatory 
behavior since they ignore spatial correlations. In higher dimensional settings it can be 
shown that mathematically more elaborate objects than one-point statistics are needed. 
However for this prototypical one-dimensional equation they suffice. A general discussion 
on the relation between Young-measure solutions and classical solutions of generalizations 



of system ([4.1|) see |28 



We summarize the results of this section. Replacing the original coordinates displacement- 
velocity {u, ii) by the new coordinates strain-strain rate {e, rj) it was possible to transform 
system (|3.1|) into the equivalent, mathematically better accessible system The 



system supports oscillating solutions and one can show that the dynamics of oscillatory 
solution can be described with time dependent Young measures which are unique solu- 
tions of ( [4.2| ). Thus the latter equation is a homogenized version of ( [4.1| ), it describes the 
macroscopic behavior of the microscopically fluctuating strain. 

5. Long-time behavior 
The discussion of the dynamics of the microkinetically regularized wave equation (|1.2|) 



shows that the stationary solutions of (|4.2|) are central for the understanding of the qual- 
itative behavior. They represent a kind of oscillatory ground state, which is to a large 
extend independent from the macroscopic dynamics. 

In the numerical simulation results in Chapter ^ one can see that the phase transfor- 
mation induces strong oscillations to the solutions. The macroscopic dynamics of these 
oscillatory solutions is very slow, this suggests the supposition that they approximate 
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stationary Young measures. This supposition is partly supported by the results of Sec- 



tion In the case p = generic solutions thermalize, i.e. {d^u, dxu){t) generates a 
nontrivial stationary Young measure as time tends to infinity. For nonzero p thermal- 
ization can be confirmed neither analytically nor numerically. The solutions still develop 
arbitrarily fine oscillations but a certain amount of macroscopic dynamics remains in the 
simulations. 

Since the microkinetic system ( |3.1| ) is Hamiltonian it might be surprising at the first 
glance, that generic solutions may converge to an equilibrium for t ^ oo. This is only 
possible by developing fine oscillations so that the macroscopically visible kinetic energy 
vanishes into uncorrelated microscopic fiuctuations. This phenomenon is denoted ther- 
malization, it is the usual behavior of solids in the absence of an exciting force. The 
fundamental laws of physics are reversible, thus the introduction of additional quantities 
like temperature, which are necessary to describe irreversible processes, only give a model 
like, statistical description of this behavior. 



5.1. Generalized stationary states. To clarify the role of system (|4.2| ) we demonstrate 
that the system is a nontrivial extension of ( [4.1|) . This is done by comparing the set of 
equilibrium points C of ( |3.1| ) with the equilibria 3^ of ( [4.2| ). It turns out that y is much 



larger than C, the latter set is not even dense in the former. The stationary points of (|3.1|) 
are characterized by the equation dx<j{dxu) = 0. Using the Neumann-boundary condition 
this equation can be integrated once and one obtains a{dxu) = 0. We denote the three 
zeros of a by Zi where e~ < zq < . The set of equilibrium points can be characterized 
as follows: For every stationary solution u* of (|3.1| ) there exists a partition (01,^0)^2) 
of Q so that dxU*{x) = Yl'i=oXni{x)zi almost everywhere. 

The equilibrium points of the extended system (^4.2|) are defined by the equation 

(5.1) {^{x, y),\/g{y) ■ (y^, -^(^(^/i) - M^)))*)^^^, = 0, 

(5.2) (p - Pdl)k = p (vl/(x, y), a(i/i))^,j,. 

for every test function g G C^(R^). By choosing g{y) = y2 we obtain the relation 

(^^J(x,y),-^iaiy^)-kix)))^^^^ = 

for every x G fi. Hence (|5.2|) implies k = 0. Thus it suffices to analyze the equation 
Vy■{<f■{Py2,-criy^)y} = 



(5.3) ^ V^(a;,?/) . J. Vi/(y) = 0, J 



1 
-1 



where H{y) = |?/2 + W{yi). Equation (|573|) expresses the fact that for fixed x G the 
function \I'(x, y) is constant on the connected components of the equipotential lines of H. 
For the convenience of the reader we give a short proof of this well known fact. 

Lemma 5.1. Let ^ , H e C\M.'^ ,R) and 

(5.4) dy,^ ■ dy,H = dy,^ ■ dy,H 
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for every y G M^. Let in addition be y^^y"^ G and G C^([0, satisfying (y9(0) = 
yW^ (^(1) = 7/(2) and f^H{ip{t)) = for every t G [0, 1]. Then ^(y^^^) = ^(y^^^) holds. 

Proof We have that ^(t/s) - ^(yi) = V^(?/(t)) ■ y{t) dt = because of (|5.4|) and since 
f^H{y{t)) = VH{y{t))-y. □ 

For double well-potentials like W{e) = — 1)^ we have to split the phase space into 
three domains in oder to represent arbitrary phase distributions: 

Pi ={ye R^\H{y) < W{zo) and y^ < Zq}, 
(5.5) P2 ={yeR'\H{y)<W{zo) and y,>zo}, 



Po =M2\(PiUP, 



2 



The sets Pi and P2 correspond approximately to the two stable phases e and e^. By 



Lemma [5.1| every stationary solution '^*{x,y) of ( |4.2|) can be expressed using three func- 
tions 

Theorem 5.2. Let \E' foe a stationary solution of ( [4.2|) which has a density function \E'* : 
f2 X ^ M>. T/ien t/iere exist three functions ai : Q x ^ i = 0, 1, 2 so t/iat 



i=0 



where Pq, Pi, P2 is the partition o/M^ defined by (|5.5|) . 

The converse of the theorem also holds true if satisfy normality constraints like 
(\E'*(x, y), l)j^gK2 = 1 for every x E Q. By taking suitable sequences {a'l)k£n we also obtain 
singular equilibria. The notation of the singular equilibrium Young measures requires the 
use of geometric measure theory, so we have decided not to characterize the full set of 
stationary solutions in this work. 

Very similar results concerning stationary Young measures for a two dimensional mag- 



neto fluid have been obtained in [0, |15|. There the most probable stationary states 
are determined by the maximum of an entropy functional. In our case the importance 
of the stationary Young measures is not clear if p > 0, cf. the results of the long-time 
integrations in Chapter 10. 



The striking difference between the stationary solutions of (|3.1|) and stationary Young 
measures is that the latter may have nonvanishing velocity components, classical equilibria 
are restricted to the line rj = 0. Since the velocity distribution of stationary Young 
measures is symmetric with respect to rj = the average velocity {"^{x, y), y2)y^^2 vanishes 
for every x E Vt. One should also note that the displacement u*{x) associated to a 
stationary Young measure \E'*(x) via 

u*{x) = I {m*{x',y),y^)y^^, dx' 







is in general not a stationary solution of (|3.1|) . This is no contradiction to the fact 
that \E'* is an equilibrium since the microscopic fluctuations are not visible from the 
macroscopic level. The discussion of the stationary solutions shows, that the continuation 

19 



of equation in the space of Young measures reveals many stationary solutions which 
have been formerly hidden. 



5.2. Thermalization. The first attempt to observe the rise of uncorrelated microscopic 
fluctuations is described in the well-known paper of E. Fermi, J. Pasta and S. Ulam 
T6| . The result, which is obtained by numerical integration of several discretized one- 



dimensional quasilinear wave equations is truly opposite of the expected behavior. It 
is found that the motion of the rod is almost periodic, especially the initial state is 
approximated up to a very high accuracy in fixed time intervals. 

We consider equation to be an example for thermalization. If p = almost every 
solution u{t,x) of (|3.1| ) converges to a stationary function u*{x) in at t — * oo, 

i.e. ( |3.1|) is a good candidate to study the transformation of macroscopic kinetic energy 
into heat. But not only u{t) converges but also {dxu{t,x),dxu{t,x)) generates a unique 
nontrivial equilibrium Young measure as time goes to infinity. Thus equation (|3.1|) 
is also an example for the dynamic creation of microstructure. The case p = is an 
important limit since the equations of motion for p > emerge from the unperturbed 
system via a compact perturbation. Thus one may hope that the essential dynamics can 
be understood by analyzing the case p = 0. 

The convergence of the solutions if p = is due to the fact that the system is completely 
integrable, hence it can be analyzed easily using canonical variables. 

Theorem 5.3. Let p = andu{t,x) be a solution of system (|3.1| ) generated by the initial 
condition {u,u)\t=o = {uq,Vq) and let the genericity assumption closuTe{G) = Q be fulfilled, 
where the subset G G Q is defined as follows: 

d_ 

dx 

Then there exists a unique nontrivial stationary Young measure \E'* G Y(i7,M^) so that 
{dxu{t,x),dxu(t,x)) generates \I'* fort — > oo. 



G = ixen 



3? {UdxVo{x)f + W{dxUo{x)) - W{e)) ^ 

em 

2' 



Using this abstract result one can easily deduce the following theorem which describes 
the long-time behavior of solutions more graphically. 



Theorem 5.4. Letu{t,x) be a solution of ( |3.1| ) and let the assumptions of Theorem \5. 5 



be satisfied. Then there exists a function u*{x) so that lim^^oo \\u{t,x) — u*{x)\\L'^(^n) + 
limt^oo x)||ioo(-f^) = 0. // in addition G = Q then there is a constant C > so that 

\\u{t,x) -U*{x)\\L^^n) + ||M(t,x)||Loo(f^) < G/t. 

For the proofs of both theorems see P7[] . 

Crucial for this result is the nonmonotonicity of a. If a is linear, the set G in Theo- 
rem 5^ is empty, the integral is constantly equal to 27ry^(3/a'. Thus the thermalization 



process is strongly related to underlying phase transformation problem. 

The discussion of the stationary solutions of the microkinetic system reveals that there 
exists a great variety of Young measure equilibria, which cannot be approximated with 
classical stationary solutions. The set of the generalized equilibria is very important for 
the dynamics, it can be shown that it attracts generic solutions in the case p = 0. The 
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simulation results show that also in the coupled case p > the solutions stay close to the 
stationary Young measures. 



6. Numerical scheme and simulation results 



Since there are no rigorous results concerning the long-time behavior if p > is consid- 
ered, we present some numerical results which should illuminate the qualitative behavior 
in the coupled case. The results indicate that the asymptotics for p > differ from the 
uncoupled case. Before we comment the individual long-time integrations we explain the 
notion of convergence of the numerical scheme. For simplicity we set /3 = 1. 

A Young-measure solution is approximated by aggregated Young measures: 



N 

^(^) = J} ^^X[hkMk+i)){^)^y^'= where y e 



k=0 i=l 



Aggregated Young measures \E''^'^ g Y'*'^(f2,M^) are constant on intervals of length h. 
For X G {hk,h{k + 1)) the probability measure \E' 



^'^[x] is a convex combination of N 
Dirac masses located at {y*'^, i = 1 . . . N}. We remark, that there are various possibilities 
of approximating nontrivial Young measures by simpler objects. For a discussion of some 
methods, see Roubicek's book ||9l. Chapter 4. 

In order to obtain useful estimates, we introduce a metric which has the property that 
I l°° 



,iYi7'^((],M2) is dense in Y(1],R2). We define 



(^2(^1,^2) 



sup (^i(x) - ^2(a;),c/) 

||Vsr||^oo(j.2,<l 



where the subscript x denotes the integration variable of the L^-norm. The metric (i2(-, ■) 
is an extension of the Wasserstein distance which metrizes vague convergence in the space 
of probability measures with first moment, see, e.g. fS^]. Of course for given \E' G Y(i7, M?) 
there is no a priori estimate for inf^h,]VgYh JV(n,R2) d2{'^^'^ , independent of \E' which can 
be used to compute in advance how h and N have to be chosen in order to obtain a 
prescribed accuracy. The propagation of this initial error can be controlled by using the 
continuity of the solution with respect to the initial value. 

Y'*'^ is invariant under the evolution defined by ( [4.2|) only if p = 0. For nonvanishing 
p the right hand side has to be projected on Y'^'^. This motivates the consideration of 
the modified evolution equation 

(6.1) vi/'^ = -V(,,,) -{m^- {ri, -(aie) - n,pB {^^ a)))'} , ^''U^o = ^0 G Y(fi, R^). 

The linear mapping tt^ is the projection onto the space of piecewise constant functions: 
'^hfi^) = -I Ixe[kh,{k+m /(^') ^ ^ [^^' + 1)^)- The dynamics defined by (|6Tl|) 



leaves Y'*'^ invariant. For initial values \E'q'^ G Y'''^(fi,M^) we obtain 



l)h)Syki(^t), 



k,i 
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where the x x 2)-vector y^^(t) solves the ordinary differential equation 



(6.2) y 



ki 



\ 



-a{yf) + 7(13 j^Y. X[kh,{k+i)h)(y{yl 

k,i 



ki\ 



This system can be integrated using an arbitrary integration scheme for ordinary differen- 
tial equations. For the error accumulated due to the modification of the spatial coupling 
an exponential bound can be derived. 

Theorem 6.1. Let (3=1 and \E'(t) and '^'^'^(t) solutions of ( |4.2| ) and ( |6.1|) satisfying 
the initial condition ^(0) = \I''^'^(0) = \I'q'^ G Y'^'^(fi,R^). Then there exist numbers 
Cii^^o) and C2{p, ||cr'||Loo(]R)) so that 

d2{^{t),^'''''{t)) < CiVp/l^gCat 

holds. 

Proof. See Appendix. □ 

To control the error, when arbitrary initial states are considered, one has to use Gronwall 
estimates for the Young-measure semiflow. This does not contribute to the understanding 



of the numerical scheme, therefore this step is omitted here, see however . 

In the numerical simulations in this work we used N = 1, i.e. only functions, not general 
Young measures. The reason we proved the convergence for the generalized scheme is, 
that we are convinced that the numerical results give a good picture of the qualitative 
dynamics although h and the time increment r are too large to see the asymptotic regime 
of the scheme. Theorem shows that it is possible to follow solutions which are strongly 
oscillating. For a proper implementation of these ideas, more sophisticated numerical 
techniques like adaptive schemes have to used, but this is beyond the scope of this work. 

6.1. Numerical results. In order to check the decay of macroscopic kinetic energy long- 
time integrations were performed. This might seem dubious since the error grows expo- 
nentially in time. At first it seems that the numerically computed solution has nothing 
to do with the exact solution for large times. The same problem occurs when integration 
schemes for ordinary differential equations are analyzed; only estimates which behave 
exponentially in time can be derived. This situation greatly improves if the viewpoint 
is slightly changed: For integration schemes for ordinary differential equations it can be 
shown that the numerically computed solution corresponds to a high precision to a slightly 
perturbed vector field, more precisely it can be shown that the error remains exponentially 
small over a time interval T = C(^) (see e.g. [0). Since we are only interested in qual- 
itative properties of the dynamics generated by certain vector fields, this weakened result 
is satisfactory. It is not our aim to give rigorous estimates for the long-time integrations, 
we just point out the ideas which are relevant if a rigorous justification of the validity of 
the results of the long-time integrations is attempted. 

As integration method for the system of ordinary differential equations the implicit 
midpoint rule with step length r = 0.1 is used. Equation (|6.2|) is integrated for different 
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values of p over the interval [0, 400], the parameter f3 is constantly kept to 1. The solutions 
which are computed in the described fashion conserve the total energy not exactly but up 
to a precision of 1% over the whole time interval. For test purposes we reduced r to the 
value 0.01, now the total energy was conserved up to an error of 0.01% and the qualitative 
results were practically unchanged. 

In order to be able to compare the results for different values of p the initial value is 
Uo{x) = \/2x + Vo{x) = in every experiment. 

The simulations have been performed with different values of h {hi = 0.005, /i2 = 0.002, 
hs = 0.001, = 0.0005). The results of runs using h = and h = h4 differ very little 
from each other, thus it can be expected that the numerical results for h = are in 
agreement with the exact behavior. 

The normahzed macroscopic kinetic energy M{t) = ^-ft'mac(^) = | /^ef^'^^^dx serves as 
an indicator for the qualitative long-time dynamics. In the figures [l^ the value M{t) is 
plotted over the time in double logarithmic scale. Time t and normalized kinetic energy 
M{t) are plotted horizontally respectively vertically. The first experiment is the uncoupled 
case p = 0. One can see clearly that M{t) converges to for t — > oo. The convergence 
rate of the computed solution is in agreement with the analytical result: M{t) < C/t 
for some constant C > 0. Since M{t) is oscillating strongly within the large integration 
interval the graph {{t,M{t)) \ t G [0,400]} resembles more to a region than a line. 

In the coupled case (p > 0) the behavior of the solutions is more complex that in the 
uncoupled case. On the one hand in the beginning the solutions develop fine structure, 
otherwise M{t) would not decay. On the other hand this process soon is hidden by chaotic 
oscillations which no longer decay. The nonvanishing of M{t) as t — cx) is due to the 
nontrivial coupling. With growing values of p the amplitude of the chaotic oscillations 
becomes larger. This can be seen in the results for p = 0.02 (Figure 0), p = 0.05 (Figure H), 
p = 0.1 (Figure ||) and p = 0.2 (Figure ^0|). These results indicate that the long-time 
behavior of the solutions in the coupled case cannot be understood by an analysis which 
is only based on considerations made for the uncoupled case. 
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Figure 9. p = 0.1, n = 2000 
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n = 2000 



7. Concluding remarks 

The high frequency oscillations due to inertial stresses are not observable at the macro 
level and represent the continuum mechanical counterpart of thermal fluctuations. Such 
stress fluctuations can overcome the stress and energy barrier (BCD) in Figure 1 when 
macroscopic stress-strain curve coincides with the Maxwell line. 

The possible scenario is as follows. After some macroscopic perturbation (e.g. pre- 
scribed velocity at the boundary), macroscopic kinetic energy transforms into microscopic 
kinetic energy, or energy of stress fluctuations. Such energy exists in each specimen before 
experiment. When under prescribed growing displacement macroscopic stress reaches the 
Maxwell stress and stress fluctuations are sufficient to overcome the barrier (BCD) in 
Figure 1, then PT will occur at constant stress. 

If for conservative systems the macroscopic kinetic energy transforms into microscopic 
kinetic energy, which is not observable macroscopically, then this can be considered as a 
macroscopic dissipation mechanism. Indeed, the numerical simulations exhibit the hys- 
teretic behavior. Such a behavior is usually interpreted in macroscopic models in terms of 
dissipative threshold which the driving force has to overcome in order to phase transfor- 
mation proceeds. The threshold value depends on the amount of the transformed phase 
as it is observed in known experiments. 

The transformation of the energetically less favorable phase into is initiated 
by a nucleation pulse. The macroscopic kinetic energy of the system is transformed 
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into microkinetic energy by the onset of strong oscillations. Tliis process decreases tlie 
macroscopic kinetic energy wfiicli is necessary to overcome tlie energy barrier between 
tlie pliases and tfie transformation process stops. To quantify this picture the notion 
of Young-measure solutions has been introduced. It has been shown that the extended 
system, which describes the dynamics of oscillatory solutions, has nontrivial equilibria 
which can store the microkinetic energy. 

We conjecture that the chaotic fluctuations are very close to a unique "most probable" 
Young-measure equilibrium in the big set of stationary Young measures. The notion of a 
probable state is meant in the sense of statistical physics where the Gibbs distributions 
are identified as maximizers of the entropy. 

Most questions concerning the dynamics of the microkinetically regularized wave equa- 
tion are still open: 

It should be possible to obtain a bound C depending on p so that the macroscopic 
kinetic energy is bounded by C{p) as time goes to infinity i.e. 

limsupKmac(i!:) <C{p). 

What happens if the more physical limit p = 1, /3 — is taken? Do the solution 
converge to solutions of the convexified equation 

where o"^ is obtained from a by replacing the nonmonotone part by the Maxwell line? This 
is only meaningful for smooth solutions. In the case of shocks the question is complicated 
by the nonuniqueness of weak solutions. 

8. Appendix 

Proof of Theorem |^. We follow the usual procedure to establish the convergence of a 
numerical scheme: Consistency and stability imply convergence. 

For the given exact Young measure solutions t ^-^ and t i-^ \E''^'^(t,x) of ([4.2|) 

and ( |6.1| ) we define x) : — : ?/ i-^ x, y) as the family of diffeomorphisms 
generated by the ordinary differential equation: 

01 (t, x) = (f)2{t, x), 02(t, x) = -((T(0i(t, x)) - pB {^{t, X, y), a(0i(t, X, y)))y^^2), 

Correspondingly the solution flow of the discretized system 

q'^=p\ / = -7r,pS(vl/^'^(t,x),a», 

is denoted by (f)^{t,x). Note, that from equation ( |4.3|) it follows immediately that \1/ and 
^'^'^ satisfy 

and 

for every test function g G C(]R^). 
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Consistency of pB with TChpB: We establish the existence of a constant C > so that 

\\{Id - 7rh)pB\\un{L^n),L'Hn)) < Ch^ 

holds. This results mainly from the smoothing property of pB. We set ujk = [hk, h{k + 
1)) C M. Firstly \\TThU — u\\L2(n) is estimated by ||92,-'u||L2(n). Using the intermediate 
value theorem for every A; G {0, . . . — 1} one obtains the existence of C,k ^ ^^k so that 
'j'xeujk ^^^^ ""^"^ ~ "^(^k) holds. Hence 



u{x') dx — u{x) 



Using this estimate we obtain 



u{^k)-uix)\< \u{s)\ds<^/\x-^k\ ( / \uis)\Ms 



\TThu{x) — u{x)\'^dx 



(8.1) 



< 



E 



A/diam(u;fc)J|M'||L2(c^^) ) dx 



1 

2 



u{x') dx' — u{x) 



< 



dx 



<h 



(8.2) <h\\u'\\L2^n)- 



The L^(n)-norm of dxpBu can be estimated by the L^(fi)-Norm of u by expanding 
u into a Fourier series and estimating the coefficients of d^pBu. Since pB is selfad- 
joint, the supremum of the modulus of the spectrum of dxpB is the operator norm in 
Lm{L^{Q),L^{Q)). Hence \\dxpB\\ < snpf,^j^^0^ = ^f. By estimate (|3), the inequal- 
ity \\{Id - vrfe)pi3||Lin(L2(Q),L2(c)) < Ch where C = ■ 



is established. Repeating the last 
step gives us the desired estimate \\{Id — iih)pB\\un{L^(n),L^(n)) < Ch"^. 



Stability: The next step is to control d2{'^o,'^{t)) using Gronwall's inequality. This 
result is a purely analytical estimate for the exact system (|4.2|) . 



My), ''^ ^sM-y2 + y2 



L2(n) 



-a o (t)s,i{,y) + pB (^0, 0- o (/)^ i) 

< !\i + iip^ii lla'iu^(M)) ■ II (^o(y), \Hy) - y|) II 

Jo 

+ / II (^o(y), 1-^(^/2) + ao + (^o(rfy), ^(^2) - fXo) I) ||i2(f,) ds + tao(l + ||pi3||) 
^0 



(8.3) 



l(^o(y), |0*(y) - y|)IL2(^) < Ce(2^+llP«lllk'lU-(M))*. 
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where C = (1 + ||pi3||)(||o-'||Loo(R) (\l/o(l/), \y\) + (Tq) and ao = |cr(0)|. Using this we can 
estimate the error due to the discretization: 



Ci2(vl>(t),v&^'^(t)) 



sup {^o,9°(l)t- 9 0(l)t) 

lip(9)<l 



< \\{^o,\4>t-^t\)\ 



L2(n) 



< 



^0, 



-a o 03^1 + a o 



)||^,(^^ + \\{Id-nh)pB\\ II (^0, o ||l2 



< / £||(v[/o,|0,- 
Jo 

+ \\nhpB\\ II (^0, W o 0^,1 - 0- o ||^2(f^) ds. 



The second term can be controlled using the stability inequality ( ^.3|) 

II (^0, o II < II (^0, W ° - ^1) II + II (^0, ^) II 

<lk'IU-w (II (^0(2/), I0t(y) - y\) II + II (^0(2/), |y|) hnn)) + (To 

<lk1U-(R)5e(2'+"''''"'l"'"^°°<«'^* + Ik'llL-(K) (^o(?/), I2/I) + ^0- 
Using Gronwall's inequality again, one gets 

^ch^e^e+\w'u^iu,\\pB\\)t . r^g(2£+iipfiii!k'ii.-(«))< + iia'iu^(M) (My), \y\) + ^ 



Setting C := C(||cr'||ioo(K)C + ||cr'||L°=(R) (^o(l/); |z/|) + <^o) yields the desired estimate 
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